Constraints on Models with a Break in the Primordial Power Spectrum 
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One of the characteristics of the "Matter Bounce" scenario, an alternative to cosmological inflation 
for producing a scale-invariant spectrum of primordial adiabatic fluctuations on large scales, is a 
break in the power spectrum at a characteristic scale, below which the spectral index changes from 
n = 1 to n = 3. We study the constraints which current cosmological data place on the location of 
such a break, and more generally on the position of the break and the slope at length scales smaller 
than the break. The observational data we use include the WMAP flve-year data set (WMAP5), 
other CMB data from BOOMERanG, CBI, VSA, and ACBAR, large-scale structure data from 
the Sloan Digital Sky Survey (SDSS, their luminous red galaxies sample). Type la Supernovae data 
(the "Union" compilation), and the Sloan Digital Sky Survey Lyman-a forest power spectrum (Lya) 
data. We employ the Markov Chain Monte Carlo method to constrain the features in the primordial 
power spectrum which are motivated by the matter bounce model. We give an upper limit on the 
length scale where the break in the spectrum occurs. 



I. INTRODUCTION 

Most parameter fitting studies assume that the spec- 
trum of primordial cosmological perturbations is a fixed 
power of the co-moving wave-number k over the range of 
scales being probed by observations. At most, a param- 
eter describing a running of the spectrum is introduced. 
However, it is possible that the primordial spectrum has 
particular features which are poorly described assuming 
a slow running superimposed over a spectrum with fixed 
slope. In fact, it was found [l[ that making an ansatz 
for the spectrum which includes special features may im- 
prove the fit to the data. 

Recently, it has been realized that the "matter bounce" 
provides a possible alternative to cosmological inflation 
for providing a spectrum of scale- invariant adiabatic fluc- 
tuations on large scales. As realized in 0, [1,|3| (see also 
0), initial perturbations starting in their vacuum state 
and exiting the Hubble radius in a matter-dominated 
phase of contraction lead to a scale-invariant spectrum 
of curvature fluctuations on super-Hubble scales in the 
contracting phase ^. If the contracting phase connects 
to an expanding phase via a non-singular bounce, then, 
as studied in [g] in the context of the Lee- Wick bounce, 
a particular realization of a non-singular bounce, then 
the co-moving curvature fluctuations pass through the 
bounce without change in the spectral slope, as long as 
scales which are large compared to the duration of the 
bounce are considered (see 0, H, @ for other studies of 
this transfer of fluctuations reaching the same conclu- 
sion). 



^ For super-Hubble scale fluctuations it is important to specify the 
gauge; we are talking about curvature fluctuations in co-moving 
gauge. 



In a contracting universe, the energy density of radi- 
ation increases faster than that of matter. Hence, if the 
initial state contained both matter and radiation, it is 
inevitable that a change in the equation of state from 
matter to radiation will occur, very much like a transfer 
from radiation to matter occurs in the expanding branch. 
As we show in the following section, scales which exit the 
Hubble radius in the radiation phase acquire an n = 3 
spectrum. If the cosmological bounce is close to time- 
symmetric, we should expect the break in the slope of 
the spectrum of primordial perturbations to occur on a 
cosmologically interesting scale, close to the scale which 
re-enters the Hubble radius during the expanding phase 
at the time teq of equal matter and radiation. 

In this Note we study the constraints which can be put 
from current observations on the transition scale where 
the spectral index changes from ri = 3 to n = 1 . Com- 
parisons of a scale dependent primordial power spectrum 
with observational data can also been found in Ref. • 

The outline of this Note is as follows: in the following 
Section we explain why in a matter bounce model a break 
in the spectrum from n = 1 on large scales to n = 3 on 
small scales will occur. In Section 3 we will discuss our 
method and the data used. Our numerical results are 
presented in Section 4, and we end with a discussion. 



II. BREAK IN THE SPECTRUM IN THE 
MATTER BOUNCE SCENARIO 

Let us review the computation of the spectrum of cos- 
mological perturbations on super- Hubble scales starting 
from vacuum initial fluctuations on sub-Hubble scales. 
We work in longitudinal gauge (see e.g. pjj for a de- 
tailed review of the theory of cosmological fluctuations 
and [l^ for a shorter overview) in which the metric is 
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given by 

ds^ = a2(ry) [(1 + 2'P)dT]^ - (1 - 2$)dx2] , (1) 

where rj is conformal time, the x denote co-moving spatial 
coordinates and <i>(x, rj) is the generahzed gravitational 
potential which carries the information about the fluctu- 
ations. As long as we are far from the bounce time, a 
good quantity to characterize the magnitude of the inho- 
mogeneities is C: the curvature fluctuation in co-moving 
coordinates which is expressed in terms of $ via 



1 



l + w 



(2) 



Ti. denoting the Hubble expansion rate in conformal time, 
a prime indicating the derivative with respect of i], and 
w — p/ p being the equation of state parameter {p and p 
are pressure and energy density, respectively). 

The variable ( is closely related to the variable v [llj 
in terms of which the action for cosmological fluctuations 
has canonical kinetic term: 



c 



(3) 



where 2(77) is a function of the cosmological background. 
If matter is modeled as a scalar field ip, then 



n 

ip' 



(4) 



If the equation of state is not changing, then z ~ a. Vac- 
uum initial conditions on sub-Hubble scales hence imply 
that on such scales 



v{k,t) = 



1 



2k 



(5) 



The equation of motion for the Fourier mode Vk of v is 



(6) 



From this it follows that v is oscillating on sub-Hubble 
scales, whereas on super-Hubble scales the evolution of v 
is determined by the function z. If the background scale 
factor evolves as 



a{t) ~ tP , 
then the super-Hubble solutions scale as 
v{r]) ~ 77" 

with 



±v 



1 1 - 3p 

2 1-p 



(7) 



(8) 



(9) 



In the matter phase we have p — 2/3 and hence 
v = —3/2 which leads to the two values for a which 
are a = 2 and a = —1. The second mode is growing 



in the contracting phase. Making use of vacuum initial 
conditions ([5]) at Hubble radius crossing, we obtain a 
scale-invariant power spectrum on super-Hubble scales 



- k^'lvkimim^mikf ^ 



const , (10) 



where the first step is the definition of the power spec- 
trum, the second step uses (jSj, the third the temporal 
scaling ([8]) of the dominant mode, and in the final step 
we have inserted the vacuum initial conditions ([5]), the 
Hubble crossing relation r}H{k) ^ k~^ and the scaling of 
the scale factor as a function of conformal time (see also 
where the calculation was done following the evolution 
of $ from initial Hubble radius crossing in the contract- 
ing phase until the post-bounce expanding phase). 

In the radiation phase (see [l^ for an earlier analysis 
of the evolution of fluctuations in a contracting radiation 
phase) we have p = 1/2 and hence — —1/2 which leads 
to the two values for a which are a — 1 and a — 0. Hence 



(11) 



where ci and C2 are constants. Thus, the dominant mode 
of V is constant on super-Hubble scales (the correspond- 
ing mode of Q is growing). A calculation analogous to 
([T0| yields 



(12) 



where we have used the constancy of v on super-Hubble 
scales and the vacuum initial conditions. This is an — 
3 spectrum. 

To summarize this section, we have shown that in- 
homogeneities originating as quantum vacuum pertur- 
bations on sub-Hubble scales and crossing the Hubble 
radius in a matter-dominated contracting phase obtain 
a scale-invariant Ug = 1 spectrum, whereas those exit- 
ing the Hubble radius during a radiation-dominated con- 
tracting phase get an = 3 spectrum. 

If our universe has emerged from a time-symmetric 
bounce, we would expect all scales which enter the Hub- 
ble radius after the time tf,q of equal matter and radi- 
ation to have a scale-invariant spectrum, whereas those 
re-entering earlier would have an = 3 spectrum. If the 
bounce is asymmetric because of the production of radi- 
ation around the bounce point, we would expect that the 
transition scale is somewhat smaller than the scale which 
re-enters the Hubble radius at teq- In the following, we 
will study what constraints which current data can give 
on the scale at which the break in the index of the power 
spectrum occurs. 
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III. METHOD AND DATA 

In our study, we perform a gfobal analysis using the 
publicly available MCMC package CosmoMC^ [11]. Our 
theory input is as follows: We assume purely adiabatic 
initial conditions. Our theory parameter space vector is: 

P = {uJb,i^c,Qs,T,ns,As,ks,ns2) , (13) 

where cob = ^bh^ and tUc = ^ch'^-, in which fib and flc 
are the physical baryon and cold dark matter densities 
relative to the critical density, is the ratio (multi- 
plied by 100) of the sound horizon to the angular di- 
ameter distance at decoupling, t is the optical depth to 
re-ionization. The remaining parameters are related to 
the primordial scalar power spectrum V^{k) which is pa- 
rameterized as 15]: 

[ lnA;(fc,o) + (n.2-l)ln(^),fc>fc,. 

(14) 

where As{kso) and A'g{kso) are defined as the amplitudes 
of initial power spectrum before and after the step, Us 
and ns2 measure the spectral index before and after the 
step in the power spectrum, a step which is located at 
the CO- moving wavenumber kg , and we match the power 
spectrum at kg during the numerical calculation in order 
to guarantee the continuity of the power spectrum. For 
the pivot scale we set kso = 0.05Mpc~^. We assume that 
the universe is spatially flat and work in the context of 
the ACDM framework with equation of state w = — 1 for 
the dark energy component. 

In our calculations, we have taken the total likeli- 
hood to be the products of the separate hkelihoods of 
CMB, SNIa, LSS and Lya. Alternatively defining = 
—2 log£, we get 

2—2 |2 ,2,2 /-ir\ 

Xtotal — XCMB + XsNIa + XlSS + XLya ■ U^J 

Using a Markov chain Monte Carlo routine the package 
then determines the best fit theory parameter values and 
the corresponding statistical error bars. The data sets 
we have used are the following: In terms of CMB data, 
we have included in our parameter fitting the WMAP5 
temperature and polarization power spectra with the rou- 
tine for computing the likelihood supplied by the WMAP 
team^. We also included some small-scale CMB measure- 
ments, such as BOOMERanG [16|], Cmjl^], VSA [l| 
and the newly released ACBAR data [l^. Besides the 
CMB information, we also combine information about 
the matter power spectrum from the "Luminous Red 
Galaxies" sample from the Sloan Digital Sky Survev[20j 



^ Available at; |http: / /cosmologist.info/cos momc / * 

3 Available at the LAMBDA website: 

|http://lambda.gsfc.nasa.gov/| 



TABLE I. Constraints on the parameters and ks from CMB, 
LSS, SNe with/without Lyman a data. We have considered 

two cases: either fixing ns2 = 3 or keeping ns2 free. The 
results are obtained by marginalizing over the other theory 
parameters. 





CMB-hLSS+SN 


CMB-fLSS-fSN-fLya 


fix ns2 — 3 


free ns2 


fix ns2 = 3 


free ns2 


ris 


057+0.0112 
u-yj' _o.oii6 


n 057+0-0115 

_o,0116 


0.960 ±0.0114 


050+0-0119 


ks 


> 0.289(2(j) 


> 0.195(2cr) 


> 1.02(2o-) 


> 1.05(2cr) 



and supernova SNIa data. In the calculation of the like- 
lihood from SNIa we have marginalized over the "nui- 
sance parameter" [2l[. The supernova data we use is 
the "Union" compilation (307 sample) [l^l ■ We also con- 
sider small scale information obtained from the Lyman-a 
forest power spectrum from the Sloan Digital Sky Survey 
(SDSS) i23rl . however we also keep its unclear systematics 
in mind |24. [25| . Furthermore, we make use of the Hub- 
ble Space Telescope (HST) measurement of the Hubble 
parameter Hq = lOOh km Mpc~^ by using a Gaus- 
sian likelihood function centered around h = 0.72 with 
standard deviation a = 0.08 [26j . 



IV. NUMERICAL RESULTS 

In this section we show the results from the global fit- 
ting to the observational data. We focus on kg , the scale 
at which the break in the power spectrum occurs, and 
ng,ng2-, the power index parameters, and we marginalize 
over the other cosmological parameters. 

The numerical results are listed in Tables I. We give 
the results obtained in two different calculations: in one 
we keep free the parameters Ug and kg and fix ng2 = 3, in 
the other we keep Ug, kg and ng2 all to be free parameters. 

In Fig. [T] we plot the one dimensional probability 
distribution of kg obtained by fitting with observational 
data. In this calculation, we have fixed the second power 
index to be ng2 — 3 and let the first index Ug as well as kg 
be free. The black solid line is obtained if we only use the 
CMB, LSS and SN la data, the red dashed line results by 
including in addition the Lyman a data. The probability 
distribution is obtained by marginalizing over the other 
cosmological parameters. 

We find that with current data we can obtain a limit 
on kg. Since the data (in particular the Lyman a data) 
tends to indicate that there is slightly less structure on 
small scales than a scale-invariant primordial power spec- 
trum would predict, and the matter bounce model pre- 
dicts more power on small scales than is obtained from a 
scale-invariant spectrum, we get a lower bound on the k- 
value of the break point. From CMB, LSS and SN data 
alone, we get kg > 0.289/i Mpc~^ at 2a C.L.. Taking 
into account the Lyman a data, the limit on kg can be 
pushed to a much smaller length scale: the 2cr limits can 
be pushed to kg > 1.02h Mpc~^. 



4 



In Fig. [2] we plot the results obtained by keeping both 
and kg free. We find there are only weak constraints 
on ns2- The data mildly prefers a red rather than a blue 
tilt (such as the "matter bounce" would predict), and 
prefers the position of the break in the power spectrum to 
be at the value ks — 19 h Mpc~^. The interpretation of 
this result is that a primordial spectrum which is almost 
scale- invariant, characterized by a fixed power index Us, 
fits the large-scale structure data well, but that there is 
less structure on scales below the break point. 




FIG. 1: Constraints on ks from the observations, assuming 
a spatially flat universe. The vertical axis is the probability, 
the horizontal axis co-moving scale. The black solid line shows 
the probability distribution for the position kst of the break in 
the spectrum obtained using the CMB -I- LSS -I- SN la data. 
The red dashed line is obtained by additionally taking into 
account the Lya data. 
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FIG. 2: Constraints on ka and the second power index ns2 
from CMB + LSS SN la Lya data. 



V. SUMMARY 

There are theoretical models such as the "matter 
bounce" which predict a transition in the power spec- 
trum from being approximately scale-invariant on length 
scales larger than some distinguished scale to being ap- 
proximately scaling as rig = 3 in shorter length scales. In 
this note, we have studied the constraints on the value 
of ks, the CO- moving momentum at which the break oc- 
curs. The current cosmological data put a 2a upper limit 
which is ks ^ o(l) hMpc^^. Future data sets will be able 
to set much tighter constraints on ks- 
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